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ABSTRACT 

Context. Activity in low-mass stars is an important ingredient in the evolution of such objects. Fundamental physical properties such 
as age, rotation, magnetic field are correlated with activity. 

Aims. We show that two components of the low-mass triple system LHS 1070 exhibit strong flaring activity. We identify the flaring 
components and obtained an improved astrometric solution for the LHS 1070 ^H/(S+C) system. 

Methods. Time-series CCD observations were used to monitor LHS 1070 in the B and Ic bands. H-band data were used to obtain 
accurate astrometry for the LHS 1070 +C) system. 

Results. We have found that two components of the triple system LHS 1070 exhibit photometric activity. We identified that compo- 
nents and S are the flaring objects. We estimate the total energy, ~ 2.0 x 10 33 ergs, and the magnetic field strength, ~ 5.5 kG, of 
the flare observed in LHS 1070 S. This event is the largest amplitude, AB > 8.2 mag, ever observed in a flare star. 
Conclusions. 
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1. Introduction 

LHS 1070 is a triple system of low -mass stars at a distance of 
7.72+0. 15 pc ( Costa et al. 2005 i. Component J[ has spectral type 
M5.5-6 and components S and C have spectral types M8.5 and 
M9-9.5, respectively (Leinert et al.|2000]). There is a great deal 



Table 1. Log of the photometric observations 



of interest in this system since the astrometric orbits (specially 
of components S+C) are well determined, making possible the 
derivation of precise masses. Seifahrt et al. (2008 ) obtained val- 
ues of 0.1 57 ±0.009 M for the combined dynamical mass of the 
S+C components and 0.272 ± 0.017 M Q for the combined dy- 
namical mass of the whole system. The masses are close to the 
H-burning limit making these objects interesting targets for de- 
tailed studies of the transition between low-mass stars and brown 
dwarfs. 

Besides the low masses, the stars in LHS 1070 present 
other interesting features: the rotational velocities are vsin/ ^ 
8, 16, 16 km s for components J[, £ and C, respectively. While 
components J[ and S present signs of activity with Ha emis- 
sion, component C does not (Re iners et al.||2007[ ). LHS 1070 
also shows intense radio emission (Berger 2006 1. Under the rea- 
sonable assumptions of coevality and spin alignm ent, |Reiners et 
al. (2007) use a Skumanich-like law (Skumanich 1972) with a 
variable breaking-law index to estimate an age of ~ 1 Gyr for 
the system. 

In this work we report on photometric evidences for strong 
activity in the J[ and S components of LHS 1070 and discuss 
the importance of these results in the context of low-mass stars. 

2. Observations 

The data were collected along an observational program on ac- 
tivity of low -mass stars that is being carried out with the facilities 



Date 


N 




Telescope 


Filter 


Jul 04, 2008 


140 


30 


1.6-m 


B 


Jul 05, 2008 


150 


30 


1.6-m 


B 


Aug 25, 2008 


653 


20 


0.6-m 


Ic 


Aug 26, 2008 


640 


20 


0.6-m 


Ic 


Aug 27, 2008 


600 


20 


0.6-m 


Ic 


Aug 28, 2008 


570 


20 


0.6-m 


Ic 


Oct 10, 2008 


100 


1 


1.6-m 


H 


Sep 01, 2009 


100 


1 


1.6-m 


H 



of Laboratorio Nacional de Astroffsica (LNA/MCT), in Brazil. 
We have selected 30 objects in the Southern Hemisphere for 
at least three differential photometry observing sessions. They 
were chosen according criteria of being low-mass objects, hav- 
ing suitable comparison stars in a small field-of-view, and be- 
ing bright enough to be observed even at the 0.6-m telescopes. 
Table Q] summarizes the characteristics of the data collected for 
LHS 1070. N is the number of individual images obtained with 
integration time t exp . The H-band images were obtained mainly 
to improve the astrometric solution for the orbital elements of 
the pair yt/(S+C). 



3. Data reduction 

The reduction of the data was done with the usual IRAF cl tasks 
and consists of subtracting a master median bias image from 
each program image, and dividing the result by a normalized 
flat-field. In the H-band, additional steps of linearization and sky 
subtraction from dithered images were used in the preparation 
of the data. We treated the photometry extraction in the B- and 
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Fig. 1. B-band differential photometry of LHS 1070 on July 04, 
2008. The light curves of LHS 1070 and of a comparison star are 
presented with full circles and stars, respectively. 

Ic/H-bands in slightly different ways. The reason for this will 
become clear in the following discussion. 

In the B-band, component J[ dominates the flux of the sys- 
tem. This allows us to extract the relevant fluxes using plain aper- 
ture photometry. In the flare event (see Figure[TJ, since additional 
light could be coming from components T> or C, we fitted a dou- 
ble 2-D Moffat function to the stellar profile. 

In the Ic-band, since in principle both & and C components 
could contribute with substantial flux, we fitted simultaneously 
three 2-D Moffat functions to the stellar profile leaving only the 
position of component J[ and the amplitudes of the three compo- 
nents as free parameters to be searched for. Essentially the same 
procedure is used for the H-band data. The distances 3[-T> and 
J[-C were obtained from the orbital elements listed in Table[3]for 
LHS 1070jg/(g+C) and th e orbital elements for LHS 1070 S/C 
from Seif ahrt et al.| ( |2008| l. Our own measurements were used 
to improve the astrometric solution for the pair Ji/CB+C) (see 
Section|4.2|). In all cases we used the amoeba routine of Press et 



al. ( 1992 1 for the fitting procedure. 

In order to increase the stability of the fits, for both bands 
we fixed the Moffat parameters derived from the profile of star 
6417-00147-1 in the Tycho catalog and used f3 = 4.765 as in 
Trujil lo et al.| ([2001'>. This star is located at Aa = 103.13 arcsec 
and AS = 157.32 arcsec from LHS 1070 M. 

Figure [T] shows the B-band differential photometry for 
LHS 1070 on July 4, 2008. As one can see, a very strong flare 
characterized by an e-folding decay time of ~ 165 sec and a fac- 
tor of ~ 80 increase in brightness with respect to quiescence was 
observed. Figure [2] shows the Ic-band light curve of an event 
with a longer time-scale (the e-folding time is ~ 1080 sec) in 
August 28th. Unfortunately, the onset of this flare was not ob- 
served. 



4. Analysis and results 

4.1. Astrometry 

Visual examination of the CCD images during the B-band 
flare suggests a significant displacement of the photocenter of 
LHS 1070 during the event. 101 pre-flare images allow the aver- 
age relative position of LHS 1070^1 with respect to a reference 
star to be measured with ~ 20 mas accuracy. As the flare pro- 
gresses, the photocenter shifts toward North and East indicating 
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Fig. 2. Ic-band differential photometry of LHS 1070 on August 
28, 2008. From top to bottom we show the aperture photometry 
for LHS 1070, the results of the 2-D Moffat function fitting for 
LHS 1070^1 and LHS 1070S+C and the aperture photometry 
for a faint comparison star in the field. 



that component J[ was not the flaring object. Figure[3]shows the 
differential positions of LHS 1 070 with respect to a reference ob- 
ject. To identify which component, S or C, was responsible for 
the event, we proceeded as follows. First, we obtained an astro- 
metric solution for the whole field using the registered pre-flare 
images. Seven stars (excluding LHS 1070) can be used for the 
astrometry. We used the IRAF ccmap task to obtain the plate 
scale (0.315 arcsec/pixel) and rotation of the images (+CP.1 with 
respect to North). The position corresponding to the flare was 
found by fitting a double 2-D Moffat function to the stellar pro- 
file in the six images closer to the flare peak. The flare position 
is shown together with the orbital solutions for LHS 1070S/C 
( Seif ahrt et al.|2008| and LHS 1070 &/(!B+C) (Table|3) in Figure 
|4j We conclude that the flaring object was component S. 

4.2. Astrometric solution and orbital fitting for LHS 1070 

Jl/(S+C) 

We used the Monte Carlo Markov Chain approach to ex- 
plore the distribution of probability of the orbital elements 
of LHS 1070.7l/(ff+ C). The data from |Seifahrt et al.| ( pOW 



Leinert et al.| ( |2001| l and our own measurements (see Table |2| 
were used in this analysis. The orbital elements distributions of 
Seifa hrt et al.| ( [2008] l were used as prior information to obtain 
the posterior distribution of the parameters. Figure [5] shows the 
marginal distributions and Table |3]presents the numerical values 
with the associated +68% uncertainties. 

It is interesting to compare the results of the orbital elements 
obtained for the LHS 1010JI/( !B+C) system with the results ob- 
tained by |Seifahrt etldj | |2008) for LHS 1070 S/C. The orbits of 
the three components are, within the uncertainties, coplanar. 

4.3. The flare in LHS 1070 £ 

As one can see in Figure [T] the B-band flare is relatively fast, 
with a rise-time of ~225s. The decay back to quiescence took 
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Table 2. Distances between the photocenters of LHS 1070 S+C 
and LHS 1070^. 



Date 


Aa (") 


AS (") 


Oct 10, 2008 
Sep 01, 2009 


0.644±0.022 
0.665±0.013 


1.344±0.022 
1.303±0.013 
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Fig. 3. The position of the photocenter of LHS 1070 with respect 
to a Tycho catalog star. The pre-fiare position was obtained us- 
ing the 101 images before the flare. Its uncertainty can be es- 
timated from the bidimensional distribution of individual posi- 
tions, and results in a positional accuracy of ~ 20 mas. The six 
images around the peak of the event define the position of the 
photocenter during the flare. The predicted positions of & and C 
components are marked with a "x" symbol considering compo- 
nent 31 at the origin. The box shown in the upper right corner 
represents the size of one pixel. 




1995 1996 1997 1998 
Epoch (yr) 



1.56 1.6 1.64 1.68 
Semimajor axis (arcsec) 




0.24 0.26 0.28 0.3 0.32 



LHS1070 I 



0.2 
Aoc(arcsec) 



-0.2 



Fig. 4. Configuration of LHS 1070 on Jul 04, 2008. The solid 
line represents the trajectory of the barycenter of LHS 1070 S+C 
around LHS 1070 which is at the origin. The dotted lines show 
the trajectories of components & and C around the barycen- 
ter of LHS 1070 S+C, assuming the two components have the 
same mass. The positions measured by |Leinert et al.| ( |2001| l and 
Seifahrt et al. ( 2008[ ) are shown as open and full squares. Our 
measurements appear as open circles. The full circle shows our 
estimated position for the flaring object. The predicted positions 
of the barycenter of LHS 1070 S+C and of components & and 
C are marked with a "x" symbol. The box shown in the upper 
right corner represents the size of one pixel of our detector 



Fig. 5. Marginal distribution of posterior probability for the or- 
bital elements of LHS 1070^/(^+0. The dashed line shows the 
median and the grey areas mark the ±68% confidence regions. 



~880 s. The duration of the event is thus ~1 100 s. The flare am- 
plitude with respect to the quiescent magnitude of component J[ 
is ~4.8 mag. However, since component S is the flaring object, 
and adopting B = 17.17 and B = 20.57 ( |Leinert et al.|2000) > as 
the quiescent B-magnitudes of components J?T and S, we obtain 
~8.2 mag as a lower limit for the flare amplitude. 
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Table 3. Orbital elements for the LHS 1070yt/(S+C) binary. 



the star ( |Hawley et al.|2003] ), we have the flux at the maximum 
of the flare given by 



Parameter 



Value 



P(yr) 
To (yr) 
fl(") 
e 

in 
no 

M LHS 1070 ( M o) 



85.1+2.7 
1996.48!°^ 
1 615+ 0029 

—0.02 1 

026 +0 021 
63 44+ - 57 
14 24 +0 - 42 
0.270±0.012 



R 2 

Fa = X-^nB A (T), 



(4) 



where T is the blackbody temperature, B,\ is the Planck func- 
tion, R t is the radius of the star and D is the distance to the 
object. Using Eq.[T]and AB ~ 4.8 for the maximum flare ampli- 
tude, the flux at the B-band pivotal wavelength is Fr ~ 7.12 x 
1CT 14 erg cirT 2 s -1 A -1 . With the aid of the mass-radius relation 
for low-mass stars IChabrier et al.|(|2009|l, and M = 0.0815 M 



4.3.1. Luminosity and total energy 

The luminosity of the flare can be estimated as follows. The 
monochromatic flux associated to the apparent quiescent mag- 
nitude B is, 



,10~(s) 



(1) 



for LHS 1070 £ ( |Leinert et al.|2000| l, the fractional area obtained 
is ~ 1.4 x 10~ 2 . The corresponding covered area and spherical 
volume are ~ 1.4 x 10 18 cm 2 and ~1.7xl0 27 cm 3 , respectively. 

Using the simplest case where the volume occupied by 
the magnetic field lines is spherical, and assuming the dimen- 
sion of the acceleration region equal to the emitting region 
( Aschwanden| [2002 ) and considering the standard model, i.e., 
the total flare energy being produced by magnetic energy decay 
(see e.g., Brown et al.|199 4; Dauphin 2007[), we can estimate the 
magnetic field strength, B, using the relation: 



where Fr = 6.32 x 10 9 erg cm 2 s 1 A 1 is the absolute flux 



corresponding to zero magnitude (Bes sell et al.|1 998 ). The qui- 
escent level observed prior to the flare is due to component 
It corresponds to B=17.17 ( |Leinert et al.|20 00). This sets a ref- 
erence level to which the flare flux may be referred to. Thus, the 
quiescent level of component J\., expressed as luminosity in the 
B-band is, 



Eboi 



If 
Sir' 



(5) 



L B = 4nD 2 F Bo AAlQ-(& ~ 5.4 x 10 27 ergs"', 



(2) 



This gives a value of B ~ 5.5 kG in the flaring region. Notice 
that this is a lower limit. |Reiners et al.| ( |2007| l obtained an av- 
erage B-field of ~ 4 kG over the entire star, from analysis of 
spectroscopic data. We would expect to have larger values of the 
magnetic field strength in a flare. 



where AA = 890 A is the B-band FWHM (Bessell 



2005) and 5. Discussion 



D — 7.72 pc is the distance to the system. The luminosity of the 
event can be obtained using Eq.|2]and the differential magnitudes 
of each point along the flare with respect to the quiescent level. 
The time-integral under the flare luminosity light curve is the 
total energy of the flare in the B-band. From this procedure we 
obtain E B ~ 9.8 X 10 31 erg. 



4.3.2. Magnetic field strength 

In order to obtain the magnetic field strength, we have to es- 
timate the bolometric energy of the event, E^. We start with 
the flare energy in the B-band, E%, and use the relation among 
optical, UV and X-ray energies obtained from multispectral ob- 
servations of flares in low-mass stars (jGershberg 2005 | l, 



20.9 E B ~ 2.0 x 10 33 ergs 



(3) 



Knowing the total energy irradiated and assuming that this 
flare has the same mechanisms as in solar flares, we can es- 
timate a lower limit for the magnetic field strength in the re- 
gion of the event. Stellar flares are caused by sudden changes 
of magnetic field strength in the stellar corona as result of the 
reconnection of magnetic field lines (see e.g., Kopp & Pneuman] 
1976 Priest|1986 1. Such changes convert magnetic potential en 



As one can see in Figure [2] component Si was responsible for 
the flare on August 28, 2008. The photometric activity detected 
in both J{ and S components is consistent with the spectroscopic 
results of |Leinert et al.| ( |2000| i. 

Photometric activity was unknown in LHS 1070 before. The 
noticeable flare in LHS 1070S shows that even close to the 
Hydrogen-burning limit, relatively old objects can still show 
significant magnetic fields. The activity level observed in the 
components of the LHS 1070 system, and the age estimated by 
Reiners et al. (2007) are in agreement with the activity lifetime- 
spectral type relation discussed by West et al. ( 2008[ ). 

In order to stress the importance of the event observed in 
LHS 1070 S, we recollect the most impulsive optical flares pre- 
viously recorded on dMe stars. Impulsive stellar flares are events 
lasting ~100-1000 s with large amplitudes. Gershberg (2005) de- 
scribes several flares among which the most impulsive was the 
one observed by de|de Jager et al.|(fl989|) on a dM5.5e star, UV 



Cet, on December 23, 1985, with amplitude of 5 mag in the 
B-band. Another impulsive flare with 7 mag amplitude at blue 
wavelengths was observed by Bond]([T976]) on a dMe star.|Stelzer 



et al.| ( 2006| ) reported a multi-wavelength observation of a flare 



ergy into plasma acceleration (see e.g., Kuperus 1 976| [Dau phin 
|2007| ). Part of the plasma is released from the star as coronal 
mass ejection and part is accelerated towards the chromosphere. 
The latter part interacts with denser plasma converting kinetic 
energy into thermal energy and is associated with optical and 
X-ray emission. According to Hawley & Fisher ( 1992 ), a black- 
body at a temperature ~ 10 4 K can be used as a raw description 
of the optical emitting region. In terms of fractional area X on 



on LP 412-31. This object has spectral type M8 and showed a 
6 mag amplitude flare in the V-band. More recently, Kowalski 
|et al. | ( 12010] ) observed a flare on the dM4.5e star YZ CMi with 
amplitude ~ 6 mag in the U-band. Thus, since the lower limit 
for the LHS 1070 S flare presented here is ~8.2 mag in the B- 
band, we conclude that this event has the largest amplitude ever 
observed in a flare star. 



Regarding the energy released in flares, Gershberg (2005) 
discusses a few events observed in BY Dra and AD Leo in the B- 
band that reached ~ 10 35 erg. However, the time-scales for these 
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events are larger than the observed by us. Besides, those objects 
have spectral types earlier than that of LHS 1070 S. This means 
that a fair comparison between the energy released in the event 
observed on LHS 1070 £ and other objects should be restricted 
to similar spectral types and similar time-scales. Using these cri- 
teria there is only one object that has a similar energy budget in 
a flare: the M8 dwarf, LP 412-31 for which E bol ~ 6 x 10 32 erg. 
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